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We present a quantitative theory for a relaxation function in a simple glass-forming model (binary 
mixture of particles with different interaction parameters). It is shown that the slowing down is 
caused by the competition between locally favored regions (clusters) which are long lived but each of 
which relaxes as a simple function of time. Without the clusters the relaxation of the background is 
simply determined by one typical length which we deduce from an elementary statistical mechanical 
argument. The total relaxation function (which depends on time in a nontrivial manner) is quanti- 
tatively determined as a weighted sum over the clusters and the background. The 'fragility' in this 
system can be understood quantitatively since it is determined by the temperature dependence of 
the number fractions of the locally favored regions. 



Two fundamental riddles in glass-forming systems 
P, Q are (i) what determines the spectacular slowing 
down of the relaxation to equilibrium when the tem- 
perature is lowered through a relatively short interval, 
and (ii) how to predict theoretically the functional forms 
of various relaxation functions. In practice one usually 
fits the data to phenomenological relaxation functions 
(e.g. the Kohlrausch- Williams- Watt KKW law) without 
any theoretical justification. Here we employ a classi- 
cal example of glass-formation, i.e. a binary mixture of 
particles with different interaction diameters, to demon- 
strate unequivocally that the slowing down is due to the 
creation of clusters of local order; these are mechani- 
cally stable and slow to relax. We present a quantitative 
computation of a (functionally non-trivial) and strongly 
temperature-dependent relaxation function by present- 
ing it as a weighted sum of cluster contributions each of 
which decays as a simple relaxation function. 
The model discussed here is the classical example Q 
of a glass-forming binary mixture of N particles in a 2- 
dimensional domain of area A, interacting via a soft 1 /r^^ 
repulsion with a 'diameter' ratio of 1.4. More or less re- 
lated models can be found in 0, @, 0, d, d, • We refer 
the reader to the extensive work done on this system 
[3, S, EH, E [is'l. The sum-up of this work is that the 
model is a bona fide glass-forming liquid meeting all the 
criteria of a glass transition. In short, the system con- 
sists of an equimolar mixture of two types of particles, 
"large" with 'diameter' (T2 = 1.4 and "small" with 'di- 
ameter' (Ti = 1, respectively, but with the same mass m. 
The three pairwise additive interactions are given by the 
purely repulsive soft-core potentials 




a,6=l,2, (1) 



where aaa = and (Jab — [o'a + o'b)/'^- The cutoff radii 
of the interaction are set at 4.5crab. The units of mass, 
length, time and temperature are m, ai, t — aiy^m/e 
and T — e/ks, respectively, with ks being Boltzmann's 
constant. The results presented below are extracted from 
molecular dynamics simulations using 30 independent 
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FIG. 1: (Color online). The relaxation function C{t) as a 
function of logarithmic time (upper panel), and of linear time 
in the lower panel. The leftmost curve (in red) pertains to 
T = 1.0 and in order to the right, the temperatures are 
T = 0.80, 0.7, 0.65, 0.6, 0.56, 054, 0.52, 0.50, 048, 0.46, 
0.44 and 0.42. Note the extreme slowing down in the range 
0.56 < T < 0.42. The black (dashed) line at T = 0.6 sepa- 
rates "simple" time dependence at higher temperature, which 
can be well fitted to a stretched exponential function from 
"non-trivial" time dependence at lower temperatures which 
is ill-fitted by any stretched exponential form. Our aim is to 
predict quantitatively the form and value of the relaxation 
function for all temperatures and times. 



systems of 4096 particles each, using the Nose-Poincare- 
Andersen thermostat [31 ■ We employ periodic boundary 
conditions on the torus. 

The relzLxation function. For the sake of this Let- 
ter we introduce a relaxation function that is made as 
follows. At time zero every large particle i in the system 
is assigned a value ci = 1 and a neighbour list consist- 
ing of its n (small or large) nearest neighbours. In time 
nearest neighbours wander into the yonder, and when the 
ith particle loses one (respectively two, three) of its near- 
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FIG. 2: (Color online). A snapshot of the system at T = 0.44. 
In colours we highlight the clusters of large particles in local 
hexagonal order. The colours have no meaning. 



est neighbours, having at that time n — 1 (respectively 
n — 2, n — 3) of the original list, it is assigned a value 
Ci{t) = 2/3, 1/3 and respectively. The relaxation func- 
tion that we monitor is C{t) = {2/N)Y,f ait). Fig. [1] 
shows how the relaxation function decays in logarithmic 
time for different temperatures as indicated in the figure 
legend. While at high temperatures (T > 0.56) the func- 
tion can be fitted to a stretched exponential, for lower 
temperatures a long tail develops (starting at T = 0.56), 
destroying any fit to a stretched exponential, as can be 
seen in the lower panel where the same data is presented 
in linear time. Note the extreme slowing down exhibited 
in the lowest eight temperatures. 

To gain a quick insight as to the nature of the slowing 
down it is advantageous to watch a movie of the decay 
of the relaxation function. Such a movie is provided as a 
supplementary information to this Letter lj| . Color cod- 
ing in green every particle whose assigned value Ci{t) > 0, 
it is obvious from the movie that a fraction of particles 
lose their green colour very rapidly, whereas the non-zero 
values of C (t) at longer times are contributed entirely by 
clusters of large particles in local hexagonal order. This 
visual observations brings us to the central theme of this 
Letter. We propose that the present glass-former be- 
longs to a generic class of glass- forming systems that have 
competing states (of crystalline order or not) which are 
locally close in energy to the ground state (which is typ- 
ically unique) . Upon rapid cooling, such systems exhibit 
patches (or clusters) of these competing states which be- 
come locally stable in the sense of having a relatively 
high local shear modulus. It is in between these clusters 



where ageing, relaxation and plasticity under strain can 
take place. A still picture of a typical realization of this 
system is provided in Fig. [2] where the locally hexagonal 
clusters of large particles are highlighted in colour (no 
meaning to different colours). For clarity of presentation 
we do not show in this figure the local hexagonal patches 
of small particles, or any other cluster of a competing 
phase. 

Cluster decomposition formula. In the rest of this 
Letter we show that this insight is the basis of a quan- 
titative theory of the relaxation function. To this aim 
we separate in our mind the clusters (here we only take 
into account the clusters of large particles in hexagonal 
local order, which are like curds) from the rest of the 
system, which we refer to as the whey. In the whey the 
relaxation process is dominated by 'defects', mainly large 
particles in pe ntag onal cages and small particles in hep- 



tagonal cages [ll|, The relaxation process necessi- 



tates a pair- wise annihilation of such defects, because of 
the Euler constraint (requiring on the torus an average of 
six neighbours per particle). If the average distance be- 
tween such 'defects' is a temperature dependent typical 
^(T), the analysis of showed that the relaxation 

of such defects is determined by a typical relaxation time 
Tw (T) which depends on the temperature according to 

r^(r)o.e''€(^)/^, (2) 

where is a constant having dimensions of fc^. Clearly, 
in order to estimate the typical distance between such 
defects we need to know how many particles belong to 
clusters. Denote the number of clusters of size s by iV^; 
then sps = sNg/N is the probability to find a large par- 
ticle in a cluster of size s. The number density of large 
particles in the whey is n^, = 1 — ^^PsS. To estimate ^ 
we denote by He and Hy^ the average enthalpy of a large 
particle inside a cluster and in the whey respectively, and 
A = He. There are w (2^ - l)/6 + 2"^/? ways to 

organize the neighbours of a large particle in the whey 
(neglecting the rare large particle in heptagonal neigh- 
bourhood), but only one way in the cluster. Thus the 
number density Uc of particles in the clusters is 

1 



This estimate, with A — 1.31 agrees very well with the 
measured value of Uc, as can be shown in Fig. [31 From 
here we can determined ^(T), the typical distance be- 
tween 'defects', as 



^-V^c (4) 

To demonstrate the usefulness of this equation we have 
computed the relaxation function C(t) only for particles 
in the whey, calling it C\o{t), and found that it can be 
excellently fitted to the simple relaxation function 



Cw(t) = aexp 



t 



(5) 
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FIG. 3: (Color online). A test of Eq. EJwith A = 1.31. In the 
inset we demonstrate Eq. ([2]) with fj, = 4.212. Here r^, was 
measured directly as the decay time of the relaxation function 
of the whey Eq. ([SJ, and S,{T) was computed from Eq. (Q. 



with a = 0.95 and (3 = 0.83. The values of t^{T) were 
extracted from such measurements at different tempera- 
tures, and plotted as \n{T^{T)) vs. ^{T)/T. The results 
are shown as the inset in Fig. [3]in very satisfactory agree- 
ment with Eq. ^ with fi = 4.212. Note that it is highly 
non-trivial that exactly the same relaxation function de- 
scribes the dynamics in the whey at all temperatures. We 
will see that this changes dramatically when the clusters 
are added to the picture. The coefficient of 0.95 repre- 
sents an almost T-independent drop in Cw{t) which re- 
sults from thermal vibrations on the scale of one r. The 
exponent 0.83 is a fit that appears to be specific to the 
present relaxation function. 

Having understood the relaxation function of the whey 
we turn now to the relaxation function of the clusters. 
Our best statistics is for relatively small clusters, and 
we demonstrate in Fig. |4] that the small clusters relax 
again as the same simple function of time as the whey, 
but with a — > a' = 0.98 and an s dependent relaxation 
time r, which can be fitted to 



Tw exp 



exp 



T 



(6) 



Note that this typical relaxation time for clusters of size 
s goes smoothly to the whey when s 0. For larger 
clusters it is more difficult to say anything precise since 
the statistics deteriorates very rapidly. Nevertheless, we 
will show now that the model ([6]) is sufficient for our 
purposes, and that we can compute the total relaxation 
function quantitatively using this simple model. To do 
so we present the relaxation function C (t) as a sum over 
the whey and the clusters in the "cluster decomposition 



FIG. 4: (Color online). Demonstration that the relaxation 
of individual clusters of size s (circles) follows the same re- 
laxation function as the whey (continuous line), but with a 
different relaxation time Ts- In blue we show s = 7atr = 0.5 
and in red s = 10 at T = 0.46. Note that we do not have 
similar fits for large clusters, and Eq. (|6]) is proposed as a 
model for all s. 



formula (CDF)" 



C{t) = nwCw{t) -I- y^p^g'exp 



rs{T) 



(7) 



It is important to note that at this point there is only 
one parameter left to fit, which is a in Eq. ([5]). The 
best fit is (J = 0.031 with which we predict the relax- 
ation functions as shown in Fig. [5] It is obvious that 
the CDF captures quantitatively both the time and the 
temperature dependence of the relaxation function. We 
submit to the reader that a comparison of data and the- 
ory of this quality had been quite rare in the subject of 
glass-forming systems. 

We strongly believe that a similar approach should be 
relevant for a whole class of glass-forming systems where 
clusters of competing phases can form upon rapid cooling 
H, [10, I2I [il, [iil . It does not mean however that 

one can automatically apply what had been done here to 
other cases. In each case the physics of the glass-former 
should be carefully understood to identify what are the 
clusters that dominate at longest times. For example in 
hydrogen bonded systems these may be compact clusters, 
fractal clusters or chains of molecules, giving difficult to 
guess formulae for as a function of the size s of the 
cluster and of the temperature T. In addition, we should 
stress that in the present case we have measured the dis- 
tribution p^; it is very desirable to derive this distribution 
from statistical mechanical first principles, as well as to 
provide a theoretical background to the fitted law 
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FIG. 5: (Color online). Comparison of the measured relax- 
ation functions with the theoretical prediction on the basis of 
the cluster decomposition formula ((TJ. From left to right the 
temperatures are T =1, 0.8, 0.65, 0.56, 0.50, 0.46 and 0.44. 
Note the quantitative agreement at all times and tempera- 
tures. The only parameter fitted here is a of Eq. ([6}. 

Notwithstanding these issues that remain for future re- 
search we propose that this example provides unequivocal 
evidence that the existence of locally favoured structure 
whose relaxation is much slower than the whey is fun- 
damental to the understanding of the phenomenology of 
the glass-transition. 
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